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SUMMARY 

A theoretical analysis -was made of the effect of aerodynamio 
hysteresis on stalling flutter. The assumption was made that the 
ahsolute magnitude of the oscillatory aerodynamic forces and 
moments eire the same at stall as at zero angle of attack hut that 
the vector magnitudes of these forces and moments are changed, this 
change being caused by the lag of aerodynamic damping and restoring 
forces behind the velocities and displacements at stall, thus 
giving rise to a hysteresis affect. Bie decrease of critical flutter 
speed at stall was thus theoretically shown. The results were 
applied to a given airfoil and correlation of the experimental and 
theoretical results was found possible by assuming that the angle 
of aerodynamic lag varies as the slope of the static-lift oiarve. 

The aerodynamio lag was shown to cause the effective torsional 
darling to decrease thereby e:q>lainlng the low values of torsional 
aerodynamio damping obtained at stall. 


IHTEODUCTIOH 

With the large increase in the use of axial-flow compressors 
and turbines, an inve'stigatlon of the vibrations of the blades of 
these units has become of Increasing importance. A significant 
part of such a study involves the type of vibration that is self- 
maintained by the continual absorption of energy from the air 
stream, namely, the flutter of turbine and oompreesor blades. 

The term flutter, as it is ordinarily used, refers to clas- 
sical flutter, a self- sustained oscillation due to the coupling of 
Inertia forces, elastic forces, damping forces, and dynamic aero- 
dynamio forces. This type of flutter usually occurs on airplane 
wings at low angles of attack when the velocity reaches a certain 
value determined by the wing design, and called the critical flutter 
speed. The theory for this type of flutter has been developed by 
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man 7 Investigators and the oscillatory aerodynainlc forces and nonents 
derived (references 1 to 3). In all these derivations, the airfoil 
is assumed to be at zero angle of attack. Close eigreement has been 
found between theory and es 2 )erlment. 

A radically diffezent tjFpe of flutter called stalling flutter, 
however, may occur on profiles at high angles of attack, such as 
blades near the stall point. This flutter occurs at much lover 
critical flutter speeds than vould be erected from classical-flutter 
theory. Two reasons exist for believing that this type of flutter 
would be more likely to occur in compressor euid turbine blades than 
classical flutter: First, compressor emd turbine blades are so 

stiff compared with ordinary airplane wings that critical speeds 
obtained by classical calculations are usually above the reasonable 
operating range, at least for subsonic airfoils. This conclusion, 
however, cannot be stated with certainty because the effect of 
oascadlhg on the oscillatory aerodynamic forces is as yet lutkoown. 
Second, the operation of compressors and turbines at high loads 
near the stall point makes the occurrence of stalling flutter more 
likely. Stalling flutter has been shown to occur on compressor 
cascades in bench tests. 

As yet, little la known about the phenomenon of stalling 
flutter. It has been observed moat frequently in connection with 
stalled propellers (reference 4). Experimental work has shown three 
possible causes (references 5 to 7). The first cause is termed 
"static Instability" (reference 5) and is associated with the fact 
that the slope of the lift curve decreases and becomes negative 
near the stall point. This change of slope can cause an instability, 
•vrtiich has been explained by Den Hartog (refearence 8) as follows: 

As the profile moves upward, the resultant motion decreases the 
effective angle of attack. If the slope of the lift curve is nega- 
tive, the lift force will Increase. As the profile moves downward, 
the effective angle of attack increases ahd the lift force decreases. 
The vibrations therefore continue to build up. This explanation in 
terms of the static lift curve is, of course, not strictly correct 
for conditions of flutter because the dynamic- lift curve will be 
different from the static- lift curve. 

The second possible cause for stalling flutter shown to exist 
(reference 5) was excitation of the airfoil by a system of K&nn£n 
vortices shed in the wake. KAmtfin vortices produce an excitation 
with frequency directly proportional to the air velocity. At 
certain discrete velocities, the Earm^ vortex frequency will 
coincide with one of the natural frequencies of the profile. When 
this coincidence of freqizencles occurs, resonant vibrations of the 
profile can be built up. 
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The third and a^rparentl^ most oonsnoii cause vas first invresti- 
gated hy Studer (reference 6) iriao carried out the first iinportant 
ez^erimental vork on the subject. Studer found that the critical 
speed of his models dropped considerably near the stall point and 
that the phase difference between the beading and torsional oscil- 
lations decreased to zero so tbat the torsicmal and bending motions 
were in phase. He also found that torsional one-degree-of -freedom 
flutter was possible, whereas it is Imown that one-degree-of -freedom 
claasioal flutter cannot occur. Studer explained the phenomenon 
occur3:dng at stall as follows: As the airfoil approaches the stall 

asagle, separation of flow occurs. As the airfoil oscillates about 
the stall point, the separation of flow is delayed until the posi- 
tion of maximum amplitude is reached. This position will be at an 
emgle of attack greater than that for the stationary airfoil. On 
the return morement, reestablishment of smooth flow is delayed to 
an angle of attack below that at which the stationary airfoil 
would stall. At the stall angle, therefore, the air forces acting 
depend on the direction of motion, ’giving rise to a hystez^sis 
loop at stall, as shown in figure 1. Studer came to the conclusion 
that this hysteresis effect, enabling energy to be absorbed from 
the air stzream, is the cause of stalling flutter. 

Studer *B work was rerifled and amplified by Victory (refer- 
ence 1 ) who used it as a basis for finding a link between stalling 
flutter and classical theory. The torsional aesrodynamic damping 
was erperimentally shown to decrease as the angle of attack 
increases and it may become negative at stall. It was then shown 
that the critical flutter speed at stall could be calculated with 
reasonable accuracy by classical theory if the usually accepted 
classical value of the torsional aerodynamic damping is replaced 
by an esperlmental value that varies with angle of attack, fre- 
quency parameter, Reynolds number, and position of the elastic 
eixis. 

A similar hysteresis phenomenon was lizvestlgated by Smllg in 
connection with aileron flutter at tremsonlc speeds. There the 
flow separation is not dtie to stall but rather to compressibility 
effects in the transonic range. The effects, however, are similar. 
The procedure used was to asstune that because of the flow separa- 
tion the danping and restoring aerodynamic forces and moments lag 
the velooltles and displacements, zrespectively. This porocedtore 
suggested a similar approach to the problem of stalling flutter. 
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The mathematical theory of the aerodyoamlo hystereala effect 
making xise of the concept of an aerodynamic angle of lag developed 
at the HACA Cleveland laboratory and the results obtained on a 
given airfoil are presented herein. 


SIMBOIfi 

The follovlng symbols are used In this analysis: 

coordinate of elastic axis measured from mldchordi 
positive towards trailing edge In units of half 
chord 


®1» ’*^1# °i# 
^2f °2 


fimctlons of reduced frequency k and angle of 
aerodynamic lag qp 


®y# ®y# aerodynamic coefficients 

b half chord used as reference unit length 

by, by, bg, bg aerodynamic coefficients 

o Hrg2 

0 i 3 stiffness of airfoil in bending per unit span 

length 

lift coefficient 

stiffness of airfoil in torsion about a per unit 
span length 

tir 

function of reduced frequency k given in ref- 
erence 1 

functions of reduced frequency k and angle of 
aerodynamic lag 

function of i*educed frequency k given in ref- 
erence 1 


Cl 

Ct 

e 

F 

G 


rH 

o> 
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I 

i 

K 

k 

L 

M 

m 

r 

s 

t 

T 

y 

yo 

a 

P 

d 


moment of Inertia alsout elastic amis per unit span 
length 

4^ 

propox*tionality constant 
reduced freq.uenoy ^ 

oscillatory aerodynamic lift force per unit spam 

oscillatory aerodynamic moment about a per unit 
span 

mass of airfoil per unit span 

mass of surrounding air cylinder per unit span, 

jtpb^ 

location of center of gravity of airfoil measured 
fzxmi a in units of half chord 

radius of gyration referred to a in units of 
half cdiord 

static mass unbalance of airfoil, mrb 
time 

critical flutter speed 

vertical displacement 

mail.TiTum bending amplitude 

static angle of attack 

phase angle between bending and torsion 

eu3gle of torsional displacement 

Tnariminn torsi<u3al amplitude 
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ratio of mass of profile per unit spen to mass of 
s;irroxindlng air cylinder per unit span, 
m o m 
“a itpl)2 

p density of svirrouQding air 

cp angle of lag between displacements or velocities 

and aerod^mamlc restoring forces or damping forces 

£0 frequency of flutter vibration 

05^ fundamental bending frequency of airfoil 

tO-t fTindamental torsional frequency of airfoil 

Dots above tbe sj^bols represent dearlvatlves vltb respect to time 


ANALYSIS 

Critical Flutter Speed and Frequency 

The oscillatory aerodynamic lift and moment are given by equa- 
tions Z7III and H of reference 1. If the exponential relations 

loot j « i(cot-^) , . . -1 1 tsb . 

y = yQe and 6 = 6q6 and the relation k » ~ atre sub- 

stituted in these equations In order to eliminate Imaginary terms, 
the equations can be written as follows: _ 

L = m^ 

M = m^b 

where 



ay’y - £^yy + (a0 -^Sy) be - Sgb 6 - ^ Ijej ^ 

byjr - £«3byy + ^bg - ^ by^ be - Wbgbe - ^ j 


( 1 ) 
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Hie second derivatire. the first derivative, arid the displace- 
ment teiais in. eq^uations (ij represent inertia terms, damping terms, 
and stiffness terms, respectively. The bending and the torsional 
displacement are, as has been already indicated, given by 


and 

e = eoei(“*“P) 

Squations (l) ‘Uien become 

L . - ffla <0® + (»e - 
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These equations represent the rector sum of the aerodynamic 
Inertia j damping, and stiffness forces and moments. 


The separation and reestahllshment of flow as the airfoil 
oscillates shout the stall point has heen shown (reference 6) to 
lag behind the actual airfoil displacement. The osolHator^ aero- 
dynamic forces and moments can he j^esolred as shown Into Inertia 
components In phase vlth the aoceleratlcms, damping ccmiponents In 
phase vlth the velocities, and stiffness components In phase with 
the displacements. If the flow separation and reestahllshment lags 
the motions, these aerodynamic force and moment components lag the 
velocities and the displacements. The essential assumption Is nov 
made that at stall the absolute values of the oscillating aero- 
dynamlo £<xcoea etnd moments are the same as at zero angle of attack 
hut because of the hysteresis «Cfeot the vector magnitudes of the 
aerodynamic restoring forces and momenta and the aerodynamic 
damping forces and momenta have so changed that the restoring forces 
and moments lag the displacements and the damping forces and 
moments lag the velocities by an eingle cp. 


Equations (2) are then modified as follows: 


L = - ]% (0^ 


®'yy'o® 


icot 




l(a)t-CF5 




-s-w 


h0o® 


I((i5t-P) 


M = - W^h 




> (3) 


+ ihQheo®^^^"^"^^ + ^ ^©0 ® 

Ic 


l(ft)t-p-cp) 


The equations 
erence l): 


of dynamic equl librium of the system are {vef- 

+ Chy + se - L = 0 
Sy + 10 + 0^0 - M » 0 


( 4 ) 
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or 


-mo? 

-so? - B.? 6o6^(“*-P> + Oteo«^(“*-P) - M - 0 


jiajfc _ Q(^ 006 ^^^"^^ - L = 0 




y (4a) 


By combinirig equations (4a) and (3) and dividing the first 
equation by a^d? and the second eq,uation by maOpb in order to 
jn a frft the coefficients dimensionless, the following eq^uations are 
obtained; •>. 


- -S. + . 0^“* + iave^^“*“‘^^ 

^■7 1% ^ mad? / ^ 


70 


.fz - S Vi(0)t-p) 




ei 


b0o = 0 




70 


> (5) 


^ mab2 mat^V 


+ lEg) 


b0Q » 0 


Let 


^Ehen 



= mcob^ 
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I = mrg2lj2 . 


mala 
S = mr'b 


^ = tXTgS = C 


s 

%"b 


= HX = e 


C-f, = ICO 4 


2 tOfc 


CO 4 




i%cA)2 ® a>2 (O' 

By BulDatituting Into eq.iiation (S) and factoring as shown 
0)^2 


+ Iji 


.0)' 


1 + 


70 + 


(ae - ^ - a) 


il + IS, 

. k 


e 


■icp 


gi((0t-p) ^ 0 


[(ty - e) + Ihy yQ + 


“v 

hfl A + o 

® k 


ru. 




0 )' 


+ ( — Z + Ihg ) e 


-icp 


gi(art-p) 130^ ^ 0 


( 6 ) 


Equations (6) ha-\^e solutions other them 0 = y = 0 if and 
only if the determinant of the coefficients of y and 6 vanishes. 




»0-| + 


c,2^.i 

CiT 


•y + - ij + iaye 

+ - [(hy - e) + ihy 

e) + lag) e-^'P 


( 7 ) 
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CD 

K 


B 7 mizltlplylng throu^ "by and substituting the Euler 

relation 

e^*^ = cos <?>+ 1 sin <p 

and hy separating real and Imaginary pairts and rearrEuaglng 

^4 / ^i\j_ v2 


(^) 




C8) 


fdiere 
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ai 

ag = 
= 

bg = 
Cl = 

Og = 




(by -e) + by i&Q - e) - bg (ay- ^) 




(bg - o)] 


OOB 


t- (sy-M’) - »y (by-©)Jain 


Cp 


914 



MA.CA EM No. E8B04 


13 


B 7 solving equation ( 8 ) for 



eu 3 d equating 





( 10 ) 


or 


(k, <P) = Fg 

For a given profile, the functions F^^ and Fg are functions 

only of the reduced frequency k and the angle of aerodynamic 
lag cp. For a given angle of aerodynamic lag cp, equation (ll) 
can he solved for reduced frequency k hy finding the intersections 
of the functions F^^ and F 2 plotted against k. The flutter 
frequency (d can then he ohtedned from equation (lO) and the 
critical flutter speed v from the relation 


(Oh 

k 


( 12 ) 


The critical flutter speeds v and the frequencies k can then he 
plotted as functions of the angle of aerodynamic lag cp. 


Aerodynamic Dancing 

The torsional aerodynamic danq)l 2 ig coefficient hg Is given hy 

ie a function of k only. After the reduced frequency k Is 
obtained as a function of the angle of aerodynarnl c lag cp hy equa- 
tion ( 11 ); the torsional aerodynamic damping hg can he calculated 
as a function of the angle of aerodynamlo lag cp. 


EESSDIffS AHP DISCIBSION 

The critical flutter speed v was calculated as a function of 
the angle of aerodynamic lag cp hy means of equations (ll) and (12) 
for wing H of reference 7, Kie constants for this wing are; 
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a a - 0,29 

'b a 3.375 Inches 

r * 0.228 

\i a 161.2 

( 0 ^ a 87.2 radleuss par second 
0 }-^ a 80.3 radians per second 

The radius of gyration -which was not given, was assimied to 

equal 0.5. The results are shown In figure 2. As the angle of 
aerodynamic lag <P Increases, the critical flutter speed v _ 
decreases. If the torsional aerodynamic damping coedfflclent 
Is now calculated by means of equation (13), the curve shown in 
figure 3 is obtained. The torsional aerodynamic damping coeffi- 
cient ’bg decreases with increasing angle of lag From fig- 
ures 2 and 3 the critical flutter speed v cai^be ob-balned as a 
function of the torsional damping coefficient b@. This relation 
is plotted in figure 4 and shows the critical flutter speed v 
deoreaslng as the torsional damping b@ decreases. The shape of 
this curve is very similar to those ob-fcalned In reference 7 for 
ratios of torsional to bending frequency close to 1. 

The results ob-balned -thvis far indicate -that the aerodynamic- 
lag affect can account for a drop in critical flutter speed v 
In the region of s-tall and also for a decrease in torsional aero- 
dynamic damping. In order to oh-baln quantitative agreement -with the 
ezperlmen-bally observed critical flutter speeds In the region of 
s-ball, however, it Is necessary -bo find a relation bet-ween the 
angle of aerodynamic lag cp eind the angle of at-taclc a. Figure 7 
of ^reference 5 shows -the -variation of (nrltloal flutter speed v 
-with angle of at-baolc a for this airfoil. By cross-plotting this 
curve -with figure 2, the emgle of aerodynamic lag 9 can be 
ob-balned as a function of angle of at-taoh a. This relation Is 
shown plotted as a solid line In figure 5. 

The question now presents Itself of whe-bher some physical 
basis can be found for the curve shown In figure 5. A relation 
be-bween angle of at-baclc a and the angle of aerodynamic lag cp 
suggests Itself -when the lift curve of the adrfoil Is considered. 

As the angle of at-baok a increases, -the sepeoatlon of flow 
increases and -bhe slope of -bhe lift curve decreases. Also, as 
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the separation of flov increases ^ the angle of aerod;naii!ic lag cp 
might he expected to increase. It would therefoze seem possihle 
that the angle of aerodynamic lag is related to the change of 
elope of the lift ovarre. The lift curve for the airfoil xised is 
given in reference 5. The assumption was nov made that the angle 
of aerodynamic lag cp is given tj 


<P 



da 


The slope of the lift curve at zero angle of attacilc is approxi- 
aiately Zn and the consteint K was calculated hy taking one 
eirhitrary point from the solid cinrve of figure 5. The value of K 
is eq^ual to 4S/x if the emgle aerodynamic lag is in 
degrees . Theref ore , 


cp = 




C|> 


- (14) 


This equation was plotted as the dashed curve of figure 5 and 
is seen to agree very veil vith the solid curve, which is based on 
erperlmental data. 

By ccanhining the dashed curve cf figure 5 with the relation 
shown in figure 2, the crltlcsal flutter speed v can nov he plotted 
as a function of aasgle of attack a. This function is shown as the 
dashed curve of figure 6, which agrees very veil with the erperi- 
mental cnarve obtained from figure 7 of ireference 5. (The slope of 
the lift curve shown in reference 5 hy means of which the dashed 
curve of fig. 5 was obtained is plotted in fig. 7 of this report.) 

The torsional aerodynamic damping h^ can nov also he plotted 
as a function of the angle o£ attack a (fig. 8). The torsional 
aerodynamic damping h^ drops sharply in the re^on of stall, the 
shape the curve being similar to the shape of the critical 
flutter-speed curve shown in figure 6. 


SIMIftET OF EBSOETS 

From the calculations of cnritical flutter speeds made for a 
glVen airfoil, it was shown that the ]^encmenon of stalling flutter 
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can at least In some cases be es^lained on the basis of an aerodynamic 
lag or hysteresis effect. The general characteristics of stalling 
flutter, namely, decreases In critical flutter speed and in effective 
aerodynamic torsional damping coefficient vith increasing angle of 
attack, were shown, to be related to the more fundamental concept of 
lagging aexodynamlc forces ai^ moments. Correlation between experi- 
mental and theoretical results was obtained by assuming that the 
aerodynamic angle of lag varies with the angle of attack In a manner 
that can he ezjilalned by the variation In the slope of the static- 
lift curve. More experimental emd theoretical work is necessary to 
more closely correlate the aerodyziamlc lag with the angle of attack. 


Flight Propulsion Eesearch Laboratoiry, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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Figure 3. - Variation of theoretical torsional aerodynamic damping coefficient 

with angle of aerodynamic lag. 
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Torsional aerodynamic damping coefficient, bg 

Figure 4. - Variation of theoretical critical flutter^ speed 
with torsional aerodynamic damping coefficient* 
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Plgiire 5 . - Variation of angle of aerodynamic lag with 

angle of attack. 
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Angle of attack, a» deg 


Pigtare 7» “ Variation of slope of lift ciirve with angle 
of attack. (Based on fig. 8, reference 5.) 


Torsional aerodynamic damping coefficient 
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Figure 8. - Variation of calculated torsional aerodynamic 
damping coefficient with angle of attack. 



